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In  this  paper  the  author  introduces  the  operator  :*  p  (y)A  with 

n 

higher  order  accuracy  for  approximation  to  the  differential  operator  D, 
where  ^  denotes  centered  difference  operator,  U  denotes  averaging 


operator , 


V»  ■  i  c.  '  -  Sm  Cm-1'  C0  -  1  • 

m»0 

a 

A  class  of  new  many-knot  spline  basis  n  :*  IP^IV))  was  suggested.  The 
smoothing  formulas 

f,  -  z  /”  0.  (“■“•)  f(  t)  dt  and  S.  f  -  I  f  A 

k,n  h  k,n  h  k,n  *•  i  k,n 
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A  NEW  WAY  FOR  CONSTRUCTING  HIGHER  ORDER 
ACCURACY  SPLINE  SMOOTHING  FORMULAS 

Dong- Xu  Qi* 

The  modern  mathematical  theory  of  spline  approximation  was  introduced  by 
I.  J.  Schoenberg  in  1946.  In  the  paper  [6}  he  studied  so-called  "B-spline 
basis".  A  B-spline  basis  can  be  normalized  in  various  ways.  One  of  them  is 
the  so  called  normalized  B-spline,  see  [2],  denoted  by  N^  k  for  the  B-spline 
function  of  degree  k  -  1  having  support  (x^,  x^+Jt) .  The  spline  smoothing 
formula  for  degree  k  -  1  to  an  arbitrary  function  f  can  be  represented  by 
S^f  »  £  f^N^  k*  This  approximation  has  been  used  in  curve  fitting 
successfully  [1],  [4]. 

In  order  to  improve  accuracy  of  the  smoothing  operator  S^,  the  author 

in  this  paper  suggests  a  new  spline  basis  denoted  ft  ,  instead  of  N.  .  . 

i,k,n 

Thus,  a  new  way  for  the  construction  of  spline  smoothing  formulas  is 

introduced.  I  prefer  calling  ^f  *  £  k  n  a  sntooth*n<T  operator  with 

grade  n  and  order  k.  In  here  when  n  “  k  Q  As  just  N^k  and  Sj^g 

* 

is  the  same  as  Sk.  since  nf  8  ^  +  <p  ,  this  is  a  class  of  many-knot 
splines. 

Concerning  higher  order  accuracy  spline  smoothing  formulas, 

I.  J.  Schoenberg  [1946]  has  already  discussed  in  [6]  and  Z.  S.  Liang  studied 
the  many-knot  spline  smoothing  [4].  My  main  attempt  in  this  paper  is  to 
suggest  a  new  way  for  constructing  them. 


* 

Department  of  Mathematics,  Jilin  University,  Changchun,  China. 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-004 1 


1*  The  smoothing  operator 


Denote  the  centered  difference  operator  by  A  ,  defined  by 

n 

\f(x>  t-  fix  +  |)  -  fix  -  |)  . 

For  simplicity  let  h  -  1,  and  A  :■  A  . 

The  B-spline  of  order  k  with  equally  spaced  knots  are  denoted  by 
and  it  can  be  represented  by 

\<x)  -  (AO*1)  6(x)  ,  (1.1) 

Hlfk(*>  #  (1.2) 

where  D-1  is  the  integral  operator/  5  is  Dirac  6-function. 

It  is  our  purpose  to  find  a  more  exact  difference  approximation  to  the 
operator  D.  I  would  like  to  choose  following  ready-made  identity. 

Fact  1.1  ([5]  p.  43) 

.  .  2.  .  2  v  ...m2  (ml)  2m+1  . . 

1+y)**'1+y  L  (-D  ( 2m-»1)  I  Y  *  (1,3) 

m*0  '  ' 

Fact  1.2  The  following  expansion 


x 


sh  x  £ 


m«0 


C  <ch 
in 


D 


m 


holds.  Set 


then 


( 2m+1 ) ! I  (2m+1 ) (2m-1 ) ...3.1  , 

.  S  n  m  |  _  1  ,\ra  ***1  n  a  « 
2m+1  m-1  '  '  ( 2m+ 1 ) 1 1 '  C0 


(1.4) 


2- 


i****^#**  nmKisi:  ■, 


Proof  Prom  (1.3) 


.  7.  .  2k  ./,  .  2  r  ,  ml  _m  2m+1 

log(y  +  /I  +  y  )  -  /i  +  y  2.  (-D  7?"T7777  2  y 


( 2m+ 1)11 


'i  ♦  y2  I  c  2” 


x  x  /  2  /  2 

Let  y  *  sh  — ,  then  ch  —  *  »1  +y  ,  x  -  2  log(y  +  *1  +  y  ).  Thus 


x  *  2  ch  -  I  Cm2”(shf) 
m«0 


-  2  ch  §  .h  f  [  c(2  .h2  f) 
m«0 


sh  x  £  C  (ch  x  -  1)m  . 

m-0  m 


Introduce  operators  E  and  UQ  defined  by 


E  f(x)  ;■  f (x  +  a)  , 


V(x)  s"  2  (f(x  +  ?>  +  f(x  "  2))f  U  *' 


U1  ' 


and  notice  the  relationships  between  those  operators  (see  [3],  p.  230) 

E  -  e°,  ch  j  -  U  , 

D  2  1  1 

,  WD  2  ~  2  2  ~  2  -r 

2  sh  —  *  e  -  e  *  E  -  D  **  A 

Use  ^  and  I  instead  of  x  and  1  in  (1.4) 
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D  ■  2  sh  r  I  C  ( ch  ~  -  I)' 

i  ra  i 

m*0 


-  I  C  <U-  I)"S 


m*0 


n 


I  C  (II  -  1)“3  +  R  , 
-  in  n 

m«0 


(1.5) 


where 


Dr  d  „  ra 

R  :■  2  sh  “  l  C  (ch  «  -  I) 
n  z  ra  2 

nr*n+1 


(1.6) 


Define  A^  as  the  first  part  of  (1.5),  i.e.. 


-r(n)  l  C  (V  -  I)”*  -  P  (|l)X  , 

a  .  n  n 

m*0 


where 


P  (U) 
n 


i  civ  - 1) 

m=*=0 


m 


i  2~j  i  (•"i>m"j<4>cm  1  <?>*2 

j-0  ra-j  3  i-0 


(1.7) 


In  the  general  case,  define 


K"'  ■■  -  vv\  • 


(1.8) 


This  is  a  collection  of  operators  approximate  to  D.  Beyond  doubt  Pn(1)  *  1, 
P0(U)  -  I. 

Fact  1.3  If  k  is  any  nonnegative  integer,  then  the  sum  of  all 

1  k 

coefficients  of  items  (Wh>  in  the  expansion  (Pn(  1^)  )  equals  to  1. 

Notice  (1.6),  the  first  term  in  Rn  for  any  h 


(1.9) 


This  implies  the  following: 


-4 


-I 


Theorem  1.1  Assume  that  f  e  C 


2n+3 


Then 


A^n>f(x)  -  Of (x)  -  2"3(n+1)Cn+1f(2n+3) 


.  t  .  .  n+1  ,  .  n+1  .  , 

where  5  e  [x - h,  x  +  -y  hj . 

Definition  We  call  the  operator  8 

grade  n  and  degree  f. 

It  is  to  be  noted  that  (A^^d  ^  is  just  as 

Here  it  is  the  smoothing  operator  of  grade  0  and 


Fact  1.4  From  ■nieorera  1 . 1 ,  if  g  @  P, 


2n+1 


on 


T<n)  -1 


n+1 


n+1 


AJ“'D  g  -  g,  all  xe  [a+yh,  b-yh), 


(C)h2(n+1> 

smoothing  operator  with 

with  X.  J.  Schoenberg's, 
degree  k» 

[a,b],  then 


2.  A  class  of  many-knot  splines 


As  has  been  already  pointed  out,  the  B-apline  My  with  equally  spaced 
knots  (h  -  1)  Is  the  result  of  the  O-th  grade  smoothing  operator  applied  to 
the  Dirac  6-function 

HJt  -  (Ad“1)MJe_1  -  (Ed-1)**  .  (2.1) 

Mow  we  use  the  smoothing  operator  1  of  grade  n  for  the 

6-function  repeatedly.  We  can  define  a  class  of  spline  functions  which  as 
more  knots  than  Ny : 

\  _  «-  <P  (»0)V  (2.2) 

K,n  n  k 

and 


‘W’  -  ■ 


k-1 


If  l  -  k,  then  fly  n(x)  -  (^<n)  which  has  knots 


_(n,k)  (n-H)k-j 

2  '  3  0,1 


2(n+1)k,  n  >  0  . 


We  often  take  t  »  k  if  without  note. 

The  following  facts  can  be  proved  easily  in  the  same  way  as  the 
corresponding  facts  for  My. 

Fact  2.1: 


(D 

(2) 

(3) 

(4) 

(5) 


Vi.00  ' 

fl.  (x)  »  0  for  all  |x|  >  , 

JCp  n  * 

Dmft.  (x)  -  (I(n>)™nv  <x),  0  <  m  <  k> 

k,n  k-ra,n 


-m 


fl  (x)  -  (A*n*  )*{xk+m”V(k+m-1)  |},  m  >  0; 


11  -  ’>  ^  0k,.lx)ax  ’  ” 

J—  * 


-6- 


-7- 


Then 


discrete  data 


3 .  Examples 


In  this  section  some  discussions  which  are  helpful  for  applications  in 
practice  will  be  given. 

From  <2.2>,  with  *  -  k.  n  -  1,  k  -  1,2, 3,4.  we  show  the  particular 
representations  as  follows: 


|x|  <  ^  , 


-  <  lx|  <  1  , 

|xl  -  1  , 


I  x  |  >  1  ; 


n2,1<X) 


50  65  .  , 

^6  36  |X|' 

42  _  49  - 

36  36  1  '' 

22  15  ,  , 

IS  +  36  |X|' 


36  36 


|x|  <  ^ 


|x|  <  1  , 

1  <  |x|  <  |  , 

|  <  |x|  <  2  , 

I  x  |  >  2  » 


°3,1(X) 


462  878  2 

432  ~  432  X  ' 

858  1584  .  ,  .  206  2 

“  T^T  1*1  +  777  x  ' 


432  432 
471  810 
432  ~  432 


,  ,  319  2 

! x |  +  432  X  ' 


627  654  ( 

432  432  1 

141  114  ... 


169  2 

432  X  ' 


|xl  <  2  ' 

\  <  |x|  <  1  , 

1  <  |x|  <  |  , 

|  <  |x|  <  2  , 


HI  .  JLii  |X|  +  x  ,  2  <  |x|  <  r  , 


ft  (x) 
4,1,x' 


-  - 


7920 

7776 

8444 

7776 

25212 

7776 

4152 

7776 

22440 

7776 

9060 

7776 

1308 

7776 

64 

7776 

0, 


20556 

2 

13059  .  , 

tttT  ,x| 

7776 

•  X 

+ 

3144 

|x| 

14268  2 

7776 

* 

7776  X 

53448 

|x| 

36036  2 

7776 

7776  x 

11328 

|x| 

7956  2 

7776 

♦ 

7776  X 

28560 

|x| 

11988  2 

7776 

7776  X 

9240 

|x| 

3132  2 

7776 

7776  X 

1128 

1 X  1 

324  2 

7776 

* 

7776  X 

48 

|v| 

12  2 

7776 

|x| 

+ 

7776  X 

9  f  JO 

7901  ,  ,3 
7776  ' 

•1--1  jxi3 
7776  ,X|  1 

J663  jx| 3 

7776  11  ' 

—  lx.3 

7776  ,X|  ' 


From  (2*2)  the  following  tables  are  givens 


Table  1: 


X 

^(x) 

0 

1 

7 

6 

+  1 

1 

1 

-  2 

2 

2 

±  1 

1 

12 

Table  2s 

°2.1 

(x) 

0 

N2(x) 

t  -  1 

1-2 

0 

1 

7 

100 

6 

72 

+  I 

1_ 

i 

35 

-  2 

2 

2 

72 

+  1 

1 

14 

12 

72 

.  3 

1 

+  — 

-  2 

72 

1*1  <  \  , 
\  <  lx|  <  1  , 

1  <  Ix|  <  |  , 

f  <  lx|  <  2  , 

2  <  |x|  <  |  , 

§  <  |x|  <  3  , 

3  <  I  x  |  <  ^  , 

|  <  Ix|  <  4  , 

|x|  >  4  . 
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Assume  four  points  in  the  plane  are  given t 


S3, 

,1f 

S4,1f 

t 

x(t) 

y(t) 

I 

X(t) 

y(t) 

1.6 

|  0.65648  j 

H  ; 

1.11870  i 

O.t 1804  | 

1.13213 

1.8 

\ 

1  0.89167 

1.08685 

0.83711 

1.09813 

2.0 

!  1*02315 

\ 

1.02315 

1.00514 

1.00514 

2.2 

\ 

;  1*08685 

j  0.89167 

!  1.09813 

|  0.83711 

2.4 

!  1.11870 

|  0.65648 

1.13213 

!  0.61804 
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